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Abstract 



Standard Model Higgs pair production at e + e _ colliders has the capability 
to determine the Higgs boson self-coupling A. I present a detailed analysis 
of the e + e~ — * ZHH and e + e~ — ► vvHH signal channels, and the relevant 
background processes, for future e + e~ linear colliders with center of mass 
energies of yfs = 0.5 TeV, 1 TeV, and 3 TeV. Special attention is given to the 
role non-resonant Feynman diagrams play, and the theoretical uncertainties 
of signal and background cross sections. I also derive quantitative sensitivity 
limits for A. I find that an e + e~ collider with y/s = 0.5 TeV can place 
meaningful bounds on A only if the Higgs boson mass is relatively close to 
its current lower limit. At an e + e _ collider with y/s = 1 TeV (3 TeV), 
A can be determined with a precision of 20 — 80% (10 — 20%) for integrated 
luminosities in the few ab _1 range and Higgs boson masses in the range mu = 
120 - 180 GeV. 
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I. INTRODUCTION 



The CERN Large Hadron Collider (LHC) is scheduled to begin operation in 2009, be- 
ginning a new era wherein the mechanism of electroweak symmetry breaking and fermion 
mass generation will be revealed and studied in detail. Although alternative mecha- 
nisms exist in theory, this is generally believed to be a light Higgs boson with mass 
114 GeV < run < 145 GeV [1-3]. More specifically, we expect a fundamental scalar sector 
which undergoes spontaneous symmetry breaking as the result of a potential which acquires 
a nonzero vacuum expectation value. The LHC will easily find a light Standard Model 
(SM) Higgs boson with moderate luminosity [4,5]. Moreover, the LHC will have the ca- 
pability to determine some of its properties [6,7], such as its fermionic and bosonic decay 
modes and couplings [8-11], including invisible decays [12] and possibly even rare decays 
to second generation fermions [13]. An e + e~ linear collider with a center of mass energy of 
350 GeV or more will be able to significantly improve these preliminary measurements, in 
some cases by an order of magnitude in precision, if an integrated luminosity of 500 fb _1 
can be achieved [14]. 

Perhaps the most important measurement after a Higgs boson discovery is of the Higgs 
potential itself, which requires measurement of the trilinear and quartic Higgs boson self- 
couplings. Only multiple Higgs boson production can probe these directly [15-17]. Several 
years ago, studies exploring the potential of the LHC, a luminosity-upgraded LHC (SLHC) 
with roughly ten times the amount of data expected in the first run, and a Very Large 
Hadron Collider (VLHC), have appeared in the literature [18-22]. There are also numerous 
quantitative sensitivity limit analyses of Higgs boson pair production in e + e~ collisions 
ranging from 500 GeV to 3 TeV center of mass energies [16,17,23-30], and for 77 — > HH [31]. 
The e + e~ studies, usually, focus on one particular Higgs mass and/or final state, only one 
center of mass energy, and, in many cases, estimate the background using a leading-log 
shower approximation. Furthermore, the effects of non-resonant Feynman diagrams are not 
taken into account. 

In this paper, I present a more thorough investigation of Higgs boson pair production 
in e + e~ collisions. I calculate the e + e~ — > ZHH and e + e~ — > vvHH [y = u e , u T ) 
signal cross sections for y/s = 0.5 TeV, 1 TeV and 3 TeV, and m H = 120 GeV, 140 GeV 
and 180 GeV, although uih = 180 GeV is disfavored by the most recent fit to electroweak 
data when direct search limits from the Tevatron experiments are taken into account [3]. 
A center of mass energy of 0.5 — 1 TeV is considered for the International Linear Collider 
(ILC) [32], whereas y/s = 3 TeV is the target energy for CERN's CLIC concept [33]. Only 
unpolarized electron and positron beams are considered. Since the cross section for the 
one-loop process e + e~ — > HH is more than one order of magnitude smaller than that 
for ZHH production [34], it is not considered here. Likewise, I will ignore the process 
e + e~ — > e + e~HH; due to the small Zee coupling, its cross section is about a factor 6 less 
than that for vvHH production [35]. I consider several final states, and estimate both 
the reducible and irreducible backgrounds using exact matrix element calculations. Special 
attention is given to the role non-resonant Feynman diagrams play, and the theoretical 
uncertainties of signal and background cross sections. All calculations are performed taking 
into account the anticipated resolution of future e + e~ detectors. Finally, I derive quantitative 
sensitivity bounds for several integrated luminosities, and compare the capabilities of the 
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considered e + e colliders with each other and those of the LHC, a luminosity upgraded 
LHC (SLHC) [18], and a Very Large Hadron Collider operating at a center of mass energy 
of 200 TeV [36]. 

The remainder of this paper is organized as follows. I first review the definition of 
the Higgs boson self-couplings and briefly discuss SM and non-SM predictions for these 
parameters in Sec. II. The methods and tools used in my calculations, together with the 
parametrized detector resolution are summarized in Sec. III. In Sees. IV and V the ZHH and 
uuHH signal channels and all relevant backgrounds are discussed. Quantitative sensitivity 
limits are calculated in Sec. VI. In Sec. VII, I finally present my conclusions. 

II. HIGGS BOSON SELF-COUPLINGS 

The trilinear and quartic Higgs boson couplings A and A are defined through the potential 

V(vh) = \m 2 H r, 2 H + \vrf H + ^Xrj%, (1) 

where t]h is the physical Higgs field, v = (v^Gf) -1 / 2 is the vacuum expectation value, and 
G F is the Fermi constant. In the SM the self-couplings are 

A = A = X SM = 7^- ( 2 ) 

Regarding the SM as an effective theory, the Higgs boson self-couplings A and A are per se 
free parameters, and S'-matrix unitarity constrains A to A < 87r/3 [37]. Since future collider 
experiments likely cannot probe A [38], I concentrate on the trilinear coupling A in the 
following. The quartic Higgs coupling does not affect the Higgs pair production processes I 
consider. 

In the SM, radiative corrections decrease A by 4—11% for 120 GeV < m# < 200 GeV [39]. 
Larger deviations are possible in scenarios beyond the SM. For example, in two Higgs doublet 
models where the lightest Higgs boson is forced to have SM like couplings to vector bosons, 
quantum corrections may increase the trilinear Higgs boson coupling by up to 100% [39]. 
In the MSSM, loop corrections modify the self-coupling of the lightest Higgs boson in the 
decoupling limit, which has SM-like couplings, by up to 8% for light stop squarks [40]. 
Anomalous Higgs boson self-couplings also appear in various other scenarios beyond the 
SM, such as models with a composite Higgs boson [41], or in Little Higgs models [42]. In 
many cases, the anomalous Higgs boson self-couplings can be parametrized in terms of higher 
dimensional operators which are induced by integrating out heavy degrees of freedom. A 
systematic analysis of Higgs boson self-couplings in a higher dimensional operator approach 
can be found in Ref. [43]. 

III. CALCULATIONAL DETAILS 

All calculations presented here have been performed at tree-level using MadEvent [44] 
which has been modified to allow for non-standard values of A. For some background calcu- 
lations it was also necessary to increase the maximum number of Feynman diagrams and/or 
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configurations allowed so that MadEvent could be used successfully. Background calculations 
involving the strong coupling constant, a s , were performed with the renormalization scale 
set to fj, — M z , where M z is the mass of the Z-boson. The SM parameters used are [45] 

G M = 1.16639 x 1(T 5 GeV" 2 , (3) 
M z = 91.188 GeV, M w = 80.419 GeV, (4) 



smH w = l-[^y a^ = y^G F si^9 w M^ (5) 



where Gf is the Fermi constant, Mw is the W mass, 6w is the weak mixing angle, and 
is the electromagnetic coupling constant in the G^ scheme. 

The basic kinematic acceptance cuts for events at the ILC and CLIC are 

E t > 15 GeV, Ej > 15 GeV, j> T > 15 GeV, 

5° < #(j,beam) < 175° , 9(j,j') > 10°, (6) 
5° < 9(£, beam) < 175° , 9(£, j) > 10°, 

where 9(k, I) is the angle between the two objects k and I, £ — e, //, E is the energy, and j) T 
is the missing transverse momentum. The f T cut is only imposed on final states where there 
are one or more neutrinos present. Likewise, the angular and energy cuts on the charged 
leptons are only imposed if there is at least one charged lepton in the event which does not 
result from the decay of a charm or bottom quark. 

In all calculations, I include minimal detector effects by Gaussian smearing of the parton 
momenta according to ILC detector expectations [46] 

M (had) = M£5 

E v ; ^E 
AE n . 0.102 

and assume that charged leptons and jets can be detected with an efficiency close to 100%. 
For the 6-tagging efficiency, e b) and the associated misidentification probabilities of light 
quark/gluon jets and c-quarks to be tagged as a 6-quarks, Pj->b an d P c ->b, I consider two 
scenarios [47,48]: 

e b = 90% with P^ b = 0.5% and P c ^ b = 10%, (8) 

and 

tb = 80% with Pj^ b = 0.1% and P c ^ b = 2%. (9) 

The energy loss in fo-jets due to 6-quark decays is taken into account via a parametrized 
function [49], and I will assume that b- and 6-quarks can be distinguished with 100% effi- 
ciency This should be a good approximation since ILC detectors are expected to be able to 
measure the electric charge of a 6-jet with an efficiency of 90% or better [48,50]. 

The identification of H — > bb and hadronic W- and Z-decays plays an important role 
in separating the ZHH and vvHH signal from the background. To identify W — > jj and 
Z — > jj decays, I require (V = W, Z) 



4 



M v -m(jj)\<8GeV, 



(10) 



where m{jj) is the di-jet invariant mass. Since the energy loss of the 6-jets distorts the 
H — > bb Breit-Wigner function and lowers the bb invariant mass, I impose 

100 GeV < m(bb) < 126 GeV for m H = 120 GeV and, (11) 
120 GeV < m(bb) < 150 GeV for m H = 140 GeV. (12) 

This captures most of the signal cross section. I also assume that, by the time Higgs pair 
production is being analyzed, the Higgs boson mass is accurately known from experiments 
at the LHC and/or an analysis of the process e + e" — > ZH. Initial state radiation and 
beamstrahlung are not included in any of the calculations presented here. 

In order to derive sensitivity limits for A, the distribution of Higgs-pair invariant mass, 
M HH , will be used. The M HH distribution is known to be sensitive to the Higgs boson 
self-coupling, in particular for small values of the Higgs-pair mass [25]. 

It is well known [16,23] that the e + e~ — > ZHH cross section decreases with increasing 
center of mass energy while a(e + e~ — > vvHH) increases. For yfs = 0.5 TeV, a{e + e~ — > 
ZHH) > a(e + e~ -> uuHH), and for ^/7s = 3 TeV, a{e + e- -> ZHH) < a{e + e- -> 
vvHH). For yfs = 1 TeV, the two cross sections are of the same order of magnitude, with 
e + e~ — > vvHH being the larger source of Higgs boson pairs. I, therefore, consider ZHH 
{yvHH) production for ^fs < 1 TeV {y/s > 1 TeV) only. For a center of mass energy of 
0.5 TeV, Higgs pair production is strongly phase space suppressed if ran > 140 GeV [21]. 
The m H = 180 GeV case, therefore, is only analyzed for y/s > 1 TeV. 

As stated before, all calculations reported here are carried out at tree level. The complete 
one-loop radiative corrections to e + e~ — > ZHH and e + e~ — > vvHH are known [51,52] to 
modify the lowest cross section by a few percent for the energy range considered in this paper. 
As I will show, for most of the final states considered in this paper, there are uncertainties 
which are significantly larger than the effect of the one-loop electroweak radiative corrections. 
Electroweak radiative corrections thus will be ignored in the following. 

IV. ZHH ANALYSIS 

The total e+e~ -> ZHH cross section at ^ = 0.5 TeV is a {ZHH) ss 0.18 fb. It 
decreases with increasing center of mass energy and Higgs boson mass. The current design 
luminosity of the ILC is [32] C = 2 x 10 34 cm _2 s _1 , corresponding to a yearly integrated 
luminosity of 200 ftr 1 . One therefore expects a total of about 40 ZHH events per year at 
such a machine, before Z- and Higgs boson decays, and experimental acceptance cuts and 
efficiencies, are taken into account. 

With an expected 6-tagging efficiency of 80% or higher, and a fairly large branching ratio 
{B{H -> 66) w 68% (33%) for m H = 120 GeV {m H = 140 GeV)), requiring HH -> bbbb 
and Z — > {£ = e, fx) offers the best chance to identify ZHH events. Unfortunately, the 
Z — > branching ratio is very small, and too few £ + £~bbbb events are left to make this 
final state viable for a measurement of the Higgs boson self-coupling. I therefore concentrate 
on the ZHH — > jjbbbb final state, where jj denotes a light jet pair consistent with originating 
from a Z-boson (i.e. satisfying Eq. (10)) and which is not tagged as a 6-pair. This explicitly 
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removes Z —> bb decays, reducing the combinatorial background and simplifying the analysis. 
Since additional jets tend to weaken the sensitivity limits for A [53], I require exactly two 
light jets, and four tagged 6-quarks in events. Two or more pairs of 6-quarks have to satisfy 
Eq. (11) or (12). 

In addition to the four Feynman diagrams contributing to the 0(a 3 ) e + e~ — > ZHH — > 
jjbbbb signal, there are approximately 8,500 0(a 6 ), 0(a 4 a 2 ) and 0(a 2 a 4 ) non-resonant and 
single Higgs resonant diagrams contributing to the same final state. Single-resonant Higgs 
and the 0(a 4 a 2 ) diagrams constitute the potentially largest background contributions. In 
Ref. [24], the role of the non-resonant and single resonant 0(a 2 a 2 ) and 0(a 4 ) diagrams 
in e + e~ — > ZHH — > Zbbbb was analyzed in the limit of a stable Z-boson. It was found 
that the contribution of the 0(a 2 a 2 ) diagrams was small compared with the signal for 
itt-h < 140 GeV. However, the background from 0(a 4 ) single Higgs production can be 
substantial, in particular for m H > 130 GeV. Misidentification of light jets and charm 
quarks also contribute to the background for jjbbbb production. 

Using the cuts and efficiencies listed in Eqs. (6) and (8), and requiring one light jet pair 
satisfying Eq. (10) and four tagged 6-jets with at least two pairs fulfilling Eq. (11), I calculate 
the cross section of the ZHH — > jjbbbb signal, the cross section for e + e~ — > jjbbbb includ- 
ing all 0(a 6 ), 0(a 4 a 2 ) and 0{a 2 a 4 ) diagrams, the e + e~ — > jjbbcc background (approxi- 
mately 7,300 0(a 6 ), 0(a 4 a 2 ) and 0(a 2 a 4 ) diagrams), the bb4j background (approximately 
15,600 0(a 6 ), 0(a 4 a 2 ) and 0(a 2 s a 4 ) diagrams), and the 0(a 4 a 2 ) e+e~ — > 6j background 
(approximately 7,600 diagrams). The signal is calculated for the SM Higgs self-coupling 
(AXhhh = 0), AXhhh = +1, and A\ HHH = -1, where 

ax h hh = ^hhh — 1 = t 1- (13) 

Asm 

The results for the Higgs pair invariant mass distribution with y/s = 0.5 TeV (1 TeV), 
and m H = 120 GeV and m H = 140 GeV are shown in Fig. 1 (Fig. 2). The solid black 
lines show the prediction for the SM e + e~ — > ZHH jjbbbb signal. Since it is impossible 
to know which 6-quark has to be paired with which 6-quark when reconstructing the Higgs 
bosons in the event, there is a combinatorial background from incorrect pairing, which is 
largest close to threshold. The dashed histogram in Fig. la shows that the combinatorial 
background indeed dominates the signal cross section close to threshold, but dies out quickly 
for larger values of Mhh- Similar results are obtained for different Higgs boson masses and 
center of mass energies. 

The dashed and dotted lines give the Mhh distribution for AXhhh = +1 and AXhhh = 
— 1, respectively. They demonstrate that the Higgs pair invariant mass distribution is sen- 
sitive to the Higgs self-coupling A, especially for small values of M H h- Positive (negative) 
values of AXhhh increase (decrease) the ZHH cross section. The magenta curves give the 
SM cross section for e + e~ — > jjbbbb including the full set of 0(a 6 ), 0(a 4 a 2 ) and 0(a 2 a 4 ) 
Feynman diagrams. For m H = 120 GeV and ^Js = 0.5 TeV and the renormalization scale 
set to n — Mz, the contribution from the single- resonant and non-resonant electroweak and 
QCD diagrams, and the interference effects between these diagrams and the ZHH signal di- 
agrams, decrease the cross section for jjbbbb production by about 10%. For ran = 140 GeV 
at the same center of mass energy, however, taking into account the single-resonant and 
non-resonant electroweak and QCD diagrams almost doubles the cross section. Since these 
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FIG. 1. The e + e~ — * ZHH — > jjbbbb cross section with four 6-tags as a function of the Higgs 
pair invariant mass, Mhh, for y/s = 0.5 TeV and a) ran = 120 GeV and b) run = 140 GeV. The 
black solid line is the prediction of the SM signal cross section. The black dashed and dotted lines 
correspond to the ZHH signal cross section for AXhhh = +1 and AXhhh = — lj respectively. 
The dashed histogram in part a) represents the combinatorial background from incorrectly assigned 
bb pairs in the signal. The magenta line shows the SM cross section for e + e~ — * jjbbbb including 
the full set of 0(a 6 ), O{o^a 2 ) and 0{a 2 s a A ) Feynman diagrams. The dashed blue line corresponds 
to the jjbbcc cross section. The cuts imposed and the efficiencies used are summarized in Eqs. (6), 
(8), and (10) - (12). 
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FIG. 2. The e + e~ — ► ZHH — ► jjbbbb cross section with four 6-tags as a function of the Higgs 
pair invariant mass, Mhh, for yfs = 1 TeV and a) mu = 120 GeV and b) ran = 140 GeV. The 
black solid line is the prediction of the SM signal cross section. The black dashed and dotted lines 
correspond to the ZHH signal cross section for AXhhh = +1 and AXhhh = — lj respectively. 
The magenta line shows the SM cross section for e + e~ — > jjbbbb including the full set of C(a 6 ), 
0{a A s a 2 ) and 0{a 2 s a A ) Feynman diagrams. The dashed blue line corresponds to the jjbbcc cross 
section. The cuts imposed and the efficiencies used are summarized in Eqs. (6), (8), and (10) 
" (12). 
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diagrams do not depend on the Higgs self-coupling, this considerably reduces the sensitivity 
to A for mu — 140 GeV. Due to the reduced phase space and the smaller H — > bb branching 
ratio, the ZHH — > jjbbbb cross section for m H = 140 GeV is about a factor 8 smaller than 
that for mil = 120 GeV. As I will show in Sec. VI, the reduction in signal cross section, 
combined with the increase of the single-resonant and non-resonant background, will make 
it very difficult to measure the Higgs self-coupling in jjbbbb production for m# > 120 GeV 
at an e + e~ collider with ^Js = 500 GeV. 

The contribution of the non-resonant diagrams to the jjbbbb cross section can, in prin- 
ciple, be reduced by imposing a tighter cut on the bb invariant mass. However, the energy 
loss of 6-quarks, combined with the finite resolution of detectors smears out the Higgs boson 
resonance over a fairly large bb invariant mass range. A more stringent m(bb) cut thus would, 
at the same time, considerably reduce the signal cross section. 

For y/s = 1 TeV and m# = 120 GeV (m# = 140 GeV), the single- resonant and non- 
resonant background diagrams lower the cross section by up to 50% (20%). The cross 
section reduction is mainly caused by the interference of resonant Z — > jj, and non-resonant 
diagrams. Note that the effect of the single- and non-resonant background can easily be 
confused with that of a non-SM Higgs self-coupling, especially for a small number of signal 
events. At \/s = 1 TeV, the reducible backgrounds are almost negligible. The SM cross 
section peaks in the region which is least sensitive to Xhhh- Given a sufficient integrated 
luminosity, this makes it possible to use the high Mhh region to normalize the cross section. 
Unfortunately, this method does not work for ^Js = 0.5 TeV, where non-standard Higgs 
self-couplings lead to a broad enhancement or reduction of the cross section, spread out over 
most of the accessible Mhh range. 

The largest reducible background to ZHH — > jjbbbb originates from jjbbcc production 
where both c-quarks are mistagged as 6-jets. This background is only significant for y/s = 
0.5 TeV and m H = 140 GeV. The bbAj background where two light jets are misidentified as 
6-jets, and the e + e~ — > 6j background where four jets are mistagged, are very small and are 
not shown in Figs. 1 and 2. 

Although the expected 6-tagging efficiency at the ILC is very high, requiring four tagged 
6-quarks reduces the observable cross section by a factor 0.66 for e& = 0.9. Since the signal 
cross section is very small, it is natural to explore whether it is advantageous to reduce the 
number of required 6-tagged jets. Reducing the number of required 6-tags from four to three 
increases the signal cross section by a factor 

i + i<i^ d4) 

ie. for e b = 0.9 by about a factor 1.44. However, the jjbbcc background increases by 
about a factor 10, while the bbAj and 6j backgrounds grow by more than two orders of 
magnitude. Furthermore, bbcjjj production now also contributes to the background. The 
Mhh distribution for ZHH — > jjbbbb with three tagged 6-quarks and \/s = 0.5 TeV (1 TeV) 
is shown in Fig. 3 (Fig. 4). To calculate the cross sections shown in these figures, I require 
one light jet pair in the mass window given in Eq. (10), one tagged bb pair and at least one bj 
combination with an invariant mass satisfying Eq. (11) or (12). For a center of mass energy 
of 0.5 TeV, the jjbbcc (dashed blue line or histogram), bbcjjj (solid blue line or histogram), 
and the bbAj background (solid red line or histogram) are all significantly larger than the 
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FIG. 3. The e + e~ — ► ZHH — > jjbbbb cross section with three 6-tags as a function of the Higgs 
pair invariant mass, Mhh, for -y/s = 0.5 TeV and a) m# = 120 GeV and b) mu = 140 GeV. The 
black solid line is the prediction of the SM signal cross section. The black dashed and dotted lines 
correspond to the ZHH signal cross section for AXhhh = +1 and AXhhh = —1, respectively. 
The magenta line shows the SM cross section for e + e~ — > jjbbbb including the full set of C(a 6 ), 
0{a^.a 2 ) and C(a^a 4 ) Feynman diagrams. The solid (dashed) blue and red lines or histograms 
correspond to the bbcjjj (jjbbcc) and bbAj cross section, respectively. The cuts imposed and the 
efficiencies used are summarized in Eqs. (6), (8), and (10) - (12). 
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FIG. 4. The e + e~ — ► ZHH — > jjbbbb cross section with three 6-tags as a function of the Higgs 
pair invariant mass, Mhh, for y^s = 1 TeV and a) mij = 120 GeV and b) m# = 140 GeV. The 
black solid line is the prediction of the SM signal cross section. The black dashed and dotted lines 
correspond to the ZHH signal cross section for AXhhh = +1 and AXhhh = — 1, respectively. 
The magenta line shows the SM cross section for e + e~ — ► jjbbbb including the full set of 0(a 6 ), 
0{a^.a 2 ) and C(q^q 4 ) Feynman diagrams. The solid (dashed) blue and red lines or histograms 
correspond to the bbcjjj (jjbbcc) and bbAj cross section, respectively. The cuts imposed and the 
efficiencies used are summarized in Eqs. (6), (8), and (10) - (12). 
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TABLE I. Cross sections in fb for e + e~ — > jjbbbb for m# = 120 GeV and m# = 140 GeV, and 
y/s = 0.5 TeV and y/s = 1 TeV. Shown are the results for the SM signal, the full set of C(a 6 ), 
0(a 4 a 2 ) and 0(a 2 a 4 ) Feynman diagrams (labeled "all"), and the total reducible backgrounds 
(labeled "bgd"), with 3, 4 and > 3 6-tags. Cross sections are listed for two sets of 6-tagging 
efficiencies and light quark/gluon jet and charm quark misidentification probabilities (see Eqs. (8) 
and (9)). 



e b = 0.9, P c ^ b = 0.1, Pj^ b = 0.005 





4 6-tags 
signal all bgd 


3 6-tags 
signal all bgd 


> 3 6-tags 
signal all bgd 


y/s = 0.5 TeV, m H = 120 GeV 
y/s = 0.5 TeV, m H = 140 GeV 


0.016 0.014 0.001 
0.002 0.004 0.002 


0.008 0.008 0.061 
0.001 0.005 0.077 


0.024 0.022 0.062 
0.003 0.009 0.079 


yfe = 1 TeV, m H = 120 GeV 
y/s = 1 TeV, m H = 140 GeV 


0.011 0.008 2- 10" 5 
0.002 0.002 2 • 10~ 5 


0.005 0.004 0.001 
0.0012 0.0015 0.0010 


0.016 0.013 0.001 
0.003 0.004 0.001 


e b = 0.8, P c ^ b = 0.02, Pj^b = 0.001 




4 6-tags 
signal all bgd 


3 6-tags 
signal all bgd 


> 3 6-tags 
signal all bgd 


y/s = 0.5 TeV, m H = 120 GeV 
v/s = 0.5 TeV, m H = 140 GeV 


0.010 0.009 3 • 10" 5 
0.001 0.002 6 • 10" 5 


0.011 0.011 0.010 
0.001 0.007 0.012 


0.021 0.020 0.010 
0.002 0.009 0.012 


y/s = l TeV, m H = 120 GeV 
y/s = l TeV, = 140 GeV 


0.007 0.005 
0.001 0.001 


0.007 0.006 2-10~ 4 
0.0017 0.0021 2 • 10~ 4 


0.014 0.011 2 • 10" 4 
0.003 0.003 2 • 10~ 4 



ZHH signal cross section. At an e + e~ collider with y/s = 1 TeV, on the other hand, these 
backgrounds are at most of the size of the signal. Thus, requiring three or more 6-tags may 
improve the sensitivity to A at a 1 TeV e + e~ collider. For y/s = 0.5 TeV, the dramatically 
increased background may well lead to a decreased sensitivity, despite the gain in the signal 
cross section. 

The results shown in Figs. 3 and 4 assume the efficiencies and misidentification probabil- 
ities of Eq. (8). A reduction of a factor 5 or more in most backgrounds can be achieved if one 
is willing to accept a slightly reduced 6-tagging efficiency; see Eq. (9). A detailed analysis of 
the sensitivity limits which one may hope to achieve requiring four tagged 6-quarks, or three 
or more tagged 6's using the efficiencies of Eq. (8) or Eq. (9) will be presented in Sec. VI. 
The cross sections for the two sets of efficiencies and misidentification probabilities, and the 
SM signal, e + e~ — > jjbbbb including the full set of 0(a 6 ), 0(aja 2 ) and 0(a 2 s a 4 ) Feynman 
diagrams, and the backgrounds are listed in Table I for three and four tagged 6-jets, together 
with the cross section for > 3 6-tags. 

One may wonder to what extent the renormalization scale uncertainty affects the cross 
sections listed in Table I. The reducible backgrounds are all of 0(a' >t s a 2 ) and thus are 
uncertain by a factor 2 ±4 or so. In the e + e~ — > jjbbbb cross section which has been calculated 
using the full set of 0(a 6 ), 0(a^a 2 ) and 0(a 2 a 4 ) Feynman diagrams, the 0{a^a 2 ) diagrams 
play an important role only for y/s = 0.5 TeV and m# = 140 GeV. In this case, I estimate 
that the scale uncertainty may change the cross section by a factor 1.5 ±1 . 

Table I shows that, for m H = 120 GeV, jjbbbb production with four tagged 6-quarks 
provides a clean, albeit low statistics, signal with a relatively small background. Including 
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the final state with three 6-tags increases the signal cross section by almost 50%. While 
the background is still small for y/s = 1 TeV, it becomes substantial for y/s = 0.5 TeV. 
The signal to background ratio for y/s = 0.5 TeV may be improved by choosing an optimal 
combination of 6-tagging efficiency and light quark/gluon jet and charm quark misidentifica- 
tion probability [28]. Furthermore, in absence of a calculation of the NLO QCD corrections 
for jjbbcc, bbcjjj and bbAj production, the background cross section is subject to a large 
renormalization scale uncertainty. For mn = 140 GeV, the signal cross section is so small 
that even with an integrated luminosity of several ab -1 , only a handful of signal events is 
produced. 

For y/s = 1 TeV, the background is already small for the efficiencies listed in Eq. (8) 
and not much is gained by choosing a different combination of 6-tagging efficiency and light 
quark/gluon jet and charm quark misidentification probability. 

The signal cross section may be further increased by relaxing the number of 6-tagged jets 
to two. In this case, one of the Higgs bosons may undergo the decay H — > cc or H — > gg. For 
m H = 120 GeV, this increases the signal cross section by about a factor 1.3. Unfortunately, 
the backgrounds increase by a much larger factor, and overwhelm the signal [54]. I therefore 
do not analyze the ZHH — > bb4j signal here. 

For van > 140 GeV, the relatively small branching ratio of B(H — > bb) m 30% makes 
it difficult to measure A in jjbbbb production, regardless of the 6-tagging efficiency. In this 
Higgs mass range, B(H — > W*W) w 50% and final states such as bbQj and ^ubbAj offer the 
possibility to more than double the number of observed ZHH events. I do not investigate 
these final states here. For mn > 2Mw, most Higgs bosons decay into a pair of W bosons. 
While there is not enough phase space for ZHH production in this region at a 500 GeV 
e + e~ collider, a small number of ZHH events may be produced for y/s = 1 TeV or above. 
However, as I will show in the following Section, the cross section for uuHH production 
is considerably larger in this region. ZHH production, therefore, is not discussed here for 
m H > 2M W . 

V. uvHH ANALYSIS 

For y/s > 1 TeV, e + e _ — > vvHH (y = u e , u^, u T ) becomes the dominant source of Higgs 
boson pairs. For v = v T only ZHH production with Z — > v^v^ v T v T contributes. For 
v — u e , a total of eight Feynman diagrams contribute, including four ZHH, Z — > u e u e and 
four vector boson fusion diagrams. In the following, I will discuss vvHH production and the 
relevant backgrounds for mn = 120 GeV, 140 GeV and mn = 180 GeV, and y/s = 1 TeV 
and 3 TeV. 

A. m n = 120 GeV 

For mn = 120 GeV, H — > bb decays dominate. I will therefore concentrate on the vvbbbb 
final state. Specifically, I require events with missing transverse momentum and four jets 
satisfying Eq. (6). A minimum of three of the jets have to be tagged as 6-jets, and there 
have to be at least two jet pairs fulfilling Eq. (11). 
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FIG. 5. The SM e + e~ — > vvHH — > cross section with three 6-tags as a function of 

the Higgs pair invariant mass, Mhh, for ^/i = 1 TeV. The black (red) solid line shows the 
ZHH — ► viDibbbb (I = fi, r) (e + e~ — > v e v e HH — > v e v e bbbb) cross section. The blue curve rep- 
resents the total e + e~ — ► vvHH — > vvbbbb cross section. The cuts imposed and the efficiencies 
used are summarized in Eqs. (6), (8), (10), and (11). 

The SM ZHH — > vi&ibbbb (I = //, r) and e + e~ — > v e v e HH v e v e bbbb Mhh differential 
cross sections for = 1 TeV are shown by the solid black and red lines in Fig. 5. The 
blue line gives the inclusive e + e~ — > vvHH — > vvbbbb cross section. For small values of 
M HH , the cross section is completely dominated by the process e + e~ — > v e v e HH, whereas 
for large Higgs pair invariant masses, ZHH production followed by Z — > ^z/; gives the 
largest contribution to the inclusive vvHH cross section. A qualitatively similar result is 
obtained at y/s = 3 TeV, however, with the e + e~ — ► v e v e HH cross section being more than 
a factor 100 larger than the ZHH rate at small values of Mhh- 

The main backgrounds to vvHH — > vvbbbb production are e + e~ — > vvbbbb at (9(afa 4 ) 
and (9(a 6 ) (about 2,300 single- and non-resonant Feynman diagrams), and vvbbcc (about 900 
0(a;?a 4 ) and C(a; 6 ) Feynman diagrams) and vvbbjj production (about 2,100 0(a 2 s a A ) and 
0(a 6 ) Feynman diagrams) where one or two charm or light quark/gluon jets are mistagged. 
Potentially dangerous backgrounds also come from bbbb and bbjj production with the missing 
transverse originating from the energy loss in 6-decays or jet energy mismeasurements. 



14 



The Higgs-pair invariant mass distribution for y/s = 1 TeV and 3 TeV for the SM 
uuHH signal (solid black curve), AXhhh = +1 (dashed black line), AXhhh — — 1 (dotted 
black line), the full set of 0(a 2 s a A ) and 0(a 6 ) uubbbb diagrams (magenta line), and the 
relevant background processes is shown in Fig. 6. In ZHH production, positive (moderately 
negative) values of AX HHH increase (decrease) the cross section. The opposite is true for 
uuHH production. For AXhhh — — 1, the cross section is strongly enhanced, while for 
AXhhh — +1 it is somewhat smaller than the SM cross section over much of the M^h 
range. Only near the HH threshold is the cross section larger than in the SM. As in the 
ZHH case, the deviations from the SM are largely concentrated at small values of M HH . 
Comparison with Fig. 2 shows that the cross section for uubbbb production is a few times 
larger than that for jjbbbb production at y/s = 1 TeV. 

When the full set of 0(a^a A ) and 0(a 6 ) diagrams instead of the uuHH signal diagrams 
is used to calculate the uubbbb cross section, the resulting differential cross section is reduced 
by 10 — 20%, except for the HH threshold region, where it is enhanced. It is thus relatively 
easy to confuse the effects of single- and non-resonant Feynman diagrams and those of 
moderately positive anomalous Higgs self-couplings. 

The uubbcc and uubbjj backgrounds are found to be at least one order of magnitude 
smaller than the signal for the efficiencies used here. Due to the renormalization scale 
uncertainty, the cross sections for these processes may vary by O(30 — 40%). The bbbb and 
bbjj backgrounds are largest for M H h close to the kinematic boundary, y/s, and drop rapidly 
for smaller values of Mhh- This is easily understood. In both cases, the invariant mass of 
the bbbb system equals the center of mass energy if the energy loss due to 6-quark decays and 
jet mismeasurements are not taken into account. While these backgrounds are substantial 
for M HH > 850 GeV at a 1 TeV machine, they decrease rapidly with increasing values of 

The bbbb and bbjj backgrounds are concentrated at large values of M^h, whereas anoma- 
lous Higgs boson self-couplings mostly affect the cross section for small values of the Higgs 
pair invariant mass. These backgrounds thus have little or no effect on the sensitivity limits 
for Xhhh- Since the region of large Higgs pair invariant masses is largely insensitive to 
Xhhh, measuring the cross section in this region would make it possible to normalize the 
cross section. Unfortunately, the e + e~ — > uuHH — > uubbbb cross section falls very quickly 
with increasing Higgs pair invariant mass, in particular for higher center of mass energies, 
resulting in a large statistical uncertainty in the large Mhh region. Accurate theoretical 
predictions of the SM uubbbb cross section thus will be indispensable for a measurement of 
the Higgs boson self-coupling in this final state. 

B. m H = 140 GeV 

The branching ratio for H — > bb drops rather quickly with increasing Higgs boson mass 
and, for m# = 140 GeV, only about 1/3 of the Higgs bosons decay into bb. At the same time, 
the H — > W*W — > 4/ branching ratio increases to about 50%. The 4/ final state consists 
of four jets with a probability of about 46% while the ivgjj (£ = e, u) final state has a 
branching ratio of about 29%. All other final states have a combined branching ratio of 25%. 
This suggests to consider the uubb4j, uubbiuejj, uu8j and uuiupGj final states in addition 
to uubbbb production. Since the decay HH — > bbAj has the largest individual branching 
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FIG. 6. The e + e~ — * vvHH — » vvbbbb cross section with three or more 6-tags as a function of 
the Higgs pair invariant mass, Mjjh-, for ran = 120 GeV and a) ^fs = 1 TeV and b) yfs = 3 TeV. 
The black solid line is the prediction of the SM signal cross section. The black dashed and dotted 
lines correspond to the vvHH signal cross section for AXhhh = +1 and AXhhh = — 1, respec- 
tively. The magenta line shows the SM cross section for e + e~ — > vvbbbb including the full set 
of O(o 6 ) and 0(a^a 4 ) Feynman diagrams. The blue and red lines correspond to the vvbbcc and 
vvbbjj cross section. The solid and dashed histograms represent the prediction of the bbbb and bbjj 
backgrounds. The cuts imposed and the efficiencies used are summarized in Eqs. (6), (8), (10), 
and (11). 
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fraction of all Higgs pair decays for m# = 140 GeV, I will consider vvbbAj production in 
addition to the vvbbbb final state in this Section. 

The results for vvbbbb production with to# = 140 GeV are shown in Fig. 7. I impose the 
same cuts as in Sec. V A, except for the cut on m(bb) which is replaced by that of Eq. (12). 
Taking the full set of 0(a 6 ) and 0(ala 4 ) Feynman diagrams into account in the calculation 
considerably enhances the cross section (magenta line) in the HH threshold region, while it 
reduces the vvbbbb rate by about 20% for larger values of Mhh- The effect of the non-signal 
diagrams contributing to vvbbbb production can be easily confused with that of a positive 
anomalous Higgs self-coupling. The vvbbcc background is substantial in the HH threshold 
region, but falls very quickly with the Higgs pair invariant mass, and becomes negligible 
for Mhh > 400 GeV. The vvbbjj background is about a factor four smaller than that 
originating from vvbbcc production for the efficiencies used here, bbjj and bbbb production 
constitute an important background for M H h > 800 GeV (M H h > 2.5 TeV) at y/s = 1 TeV 
(v/i = 3TeV). 

When calculating the cross section for vvHH — > vvbbAj production care has to be 
taken. In addition to H — > W*W — > Aj, H — > Z*Z — > Aj also contributes, albeit with 
a much smaller branching ratio (B(H — > Z*Z) m 10% vs. B(H — > W*W) « 50% for 
itlh = 140 GeV). To select vvHH — > vvbbAj events, I require, in addition to the standard jet 
and missing transverse momentum cuts of Eqs. (6), two tagged 6-quarks satisfying Eq. (11) 
and four non-tagged jets with an invariant mass 

\m H - rn(Aj)\ < 20 GeV. (15) 

One un-tagged jet pair has to be consistent with originating from a W decay (see Eq. (10)). 

The main reducible backgrounds originate from vv6j production where two jets are 
mistagged as 6-quarks, vvccAj production where both charm quarks are misidentified as 
&'s, and from e + e _ — > bbAj with the missing transverse momentum originating from jet 
mismeasurements and the energy loss arising from 6-decays. Furthermore, non-resonant 
vvbbAj production constitutes a potentially dangerous irreducible background. 

With the exception of e + e~ — > bbAj, these are processes with eight particles in the final 
state. I have attempted to calculate the e + e~ — > vvbbAj cross section including the full 
set of 0(aga A ), O(a 2 s o. & ) and 0(a 8 ) diagrams using several publically available programs. 
Generating the contributing Feynman diagrams via MadGraph [44] resulted in well over 10 5 
Feynman diagrams and took more than 200 hours of CPU time on a 3ghz Xeon workstation. 
The evaluation of the cross section for this process using MadEvent would require comput- 
ing resources significantly larger than those available, and therefore was not attempted. 
Like MadEvent, the current version of Sherpa [55] is based on a Feynman diagrammatic 
approach [56], and thus is expected to be too slow to calculate the cross section of 2 — > 8 
processes. In the future, Sherpa will use [57] COMIX [58], a new matrix element generator 
which is based on color dressed Berends-Giele recursive relations [59]. This should allow for 
a much faster evaluation of matrix elements, and thus make it possible to calculate the cross 
sections of processes with > 8 particles in the final state. 

For processes with many particles in the final state, WHIZARD [60] and HELAC-PHEGAS [61] 
are potential alternatives to MadEvent and Sherpa. WHIZARD uses 0'Mega [62] which im- 
plements an algorithm that collects all common sub-expressions in the sum over Feynman 
diagrams contributing to a given scattering amplitude at tree level. HELAC-PHEGAS calculates 
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FIG. 7. The e + e~ — * vvHH — » uvbbbb cross section with three or more 6-tags as a function of 
the Higgs pair invariant mass, Mjjh-, for ran = 140 GeV and a) ^fs = 1 TeV and b) yfs = 3 TeV. 
The black solid line is the prediction of the SM signal cross section. The black dashed and dotted 
lines correspond to the vvHH signal cross section for AXhhh = +1 and AXhhh = — 1, respec- 
tively. The magenta line shows the SM cross section for e + e~ — > vvbbbb including the full set 
of 0(a 6 ) and 0(a^a 4 ) Feynman diagrams. The blue and red lines/histograms correspond to the 
uubbcc and vvbbjj cross section. The solid and dashed histograms represent the prediction of the 
bbbb and bbjj backgrounds. The cuts imposed and the efficiencies used are summarized in Eqs. (6), 
(8), (10), and (12). 



18 



matrix elements using recursive Schwinger- Dyson equations. I have attempted to compute 
the e + e~ — > vvbbAj cross section using WHIZARD and HELAC-PHEGAS. The compilation of the 
WHIZARD code was terminated without a result after more than 40 hours of CPU time on 
a 3ghz Xeon workstation. The HELAC-PHEGAS code compiled successfully, but failed to run 
with an error in an underlying basic linux library. Another program which may be able to 
handle a process such as e + e _ — > vvbbAj is carlomat [63], which, however, is not publically 
available yet. 

One can argue that a substantial portion of the contribution of the non-resonant 0(a 2 a 6 ) 
and 0(a 8 ) diagrams 1 to the vvbbAj cross section originates from the off-shell W* — > jj jet 
pair. This suggests that calculating the e + e~ — > vvWjjbb cross section with W — > jj, 
or the vvHAj cross section with H — > bb, may be sufficient for an estimate of how non- 
resonant diagrams affect the vvbbAj rate. Results based on a calculation of the process 
e + e~ — > vvWjjbb with W — > jj (about 7,000 0(a 2 a 5 ) and 0(a 7 ) diagrams) are shown 
in Fig. 8. The black solid line is the prediction of the SM signal cross section; the black 
dashed and dotted lines correspond to the vvHH signal cross section for A\ HHH = +1 
and AXhhh — — 1, respectively. With the acceptance cuts imposed here, the vvbbAj cross 
section is approximately as large as that for vvbbbb production. The e + e _ — > vvWjjbb, W — > 
jj, cross section is shown by the red line. While the non- resonant QCD and electroweak 
diagrams included in e + e~ — > vvWjjbb, W — > jj, have a fairly small effect at small values of 
M HH , they are seen to reduce the differential cross section by about 30% in the large Higgs 
pair invariant mass region. The non-resonant diagrams not included in vvWjjbb, W — > jj 
may well affect the cross section to a similar degree. 

I have not attempted to compute the vvccAj and vv&j backgrounds, suspecting that one 
would run into problems similar to those encountered for e + e~ — > uubbAj. I expect that 
these backgrounds are small compared with the signal, as in the case of the analogous uubbcc 
and vvAj backgrounds in vvbbbb production. 

The e + e^ -> bUj background (with about 15,000 0(a 6 ), C(a 2 s a 4 ) and 0(a 4 s a 2 ) Feyn- 
man diagrams contributing) with the missing transverse momentum originating from jet 
mismeasurements and the energy loss arising from 6-decays was found to be very small for 
the cuts imposed. 

C. m H = 180 GeV 

For run = 180 GeV, almost all Higgs bosons decay into a pair of W^-bosons {B(H — > 
W + W~) rs 93%). Subsequent W decay then leads to HH — > t^vfij (I = e, fj) with a 
branching ratio of about 24%, or HH — > 8j with a branching fraction of ~ 19%. Since 
the individual branching ratios for all other final states are significantly smaller, I shall 
concentrate on HH — > t^vSj and HH — > 8j here. 

If one of the four lU-bosons decays leptonically, the final state consists of one charged 
lepton, six jets and missing transverse momentum which originates from the three neutrinos 
in the event. The main backgrounds originate from single resonant and non- resonant e + e~ — > 



There are no C(afa 4 ) diagrams contributing to the process. 
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FIG. 8. The e + e~ — ► vvHH — > vvbbAj cross section as a function of the Higgs pair invariant 
mass, Mhh, for m# = 140 GeV and a) \/i = 1 TeV and b) = 3 TeV. The black solid line is 
the prediction of the SM signal cross section. The black dashed and dotted lines correspond to the 
vuHH signal cross section for AXhhh = +1 and AXhhh = — 1, respectively. The red line shows 
the SM cross section for e + e~ — ► vDWjjbb, W — ► jj, including the full set of 0(a 7 ) and 0(a^a 5 ) 
Feynman diagrams. The cuts imposed and the efficiencies used are summarized in Eqs. (6), (8), 
(10), (11), and (15). 
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vv^ufij diagrams and W ± 6j production (with approximately 21,000 0{a z af) Feynman 
diagrams contributing). In order to identify vvHH events in the £ ± ^ T 6j final state, I 
require, in addition to the standard lepton, jet and f T identification cuts of Eq. (6), three 
or more jet pairs which satisfy Eq. (10) with two of the jet pairs having an invariant mass 
in the range 

160 GeV < m([jj][jj}) < 200 GeV. (16) 

Since the invariant mass of the Higgs pair cannot easily be reconstructed, I consider the 
distribution of the invariant mass of the six jets, Mqj, instead. 

Although a full calculation of the 2 — > 10 process e + e~ — > vv^vSj is currently not 
feasible, it is still possible to get an idea of how strongly non-resonant diagrams affect the 
6j invariant mass distribution by calculating the cross sections of the processes e + e~ — > 
uuiuijjH with H — > W + W~ — > 4j (with approximately 1,300 0(a 7 ) Feynman diagrams 
contributing), and e + e~ — > vvAjH with H — > W + W~ — » ivejj (with approximately 20,000 
0(a 7 ) and 0{a 5 a 2 s ) Feynman diagrams contributing). In Fig. 9, I show the SM e + e _ — > 
u e u e HH — > VeU^vfij M 6 j differential cross section for the signal at y/s = 1 TeV, together 
with the results for h> e v e ^i>ijjH , H — > W + W~ — > 4j production (red line), and e + e _ — > 
v e is e 4jH, H — > W + W~ — > (blue curve). Similar results are obtained for = 3 TeV, 

and for e + e~ — > UiUiHH — > ViUi^vfij with / = /i, r. 

The non-resonant diagrams in e + e~ — > u e i? e £ ± i'£jjH, H — > W + W~ — > 4j are found to 
significantly enhance the v e v e HH cross section near threshold. This is to be expected since 
the £z/£ invariant mass cannot be constrained; the presence of three neutrinos in the final 
state makes it impossible to use the tj> T transverse mass to reduce the background from 
non-resonant diagrams. The non-resonant diagrams in u e u e 4jH, H W + W~ — > ivgjj 
production, on the other hand, reduce the cross section by a factor 1.5 — 2 for the cuts 
imposed. An estimate of the effect of the full set of non-resonant diagrams may be obtained 
by averaging the z/ e z/ e £ ± z/^jjif, H — * W + W~ — > 4j and v e i? e 4jH, H — > W + W~ — > 
cross sections which is shown by the magenta line in Fig. 9. 

The estimate obviously ignores a large number of non-resonant Feynman diagrams. As 
demonstrated by the blue line in Fig. 9, non-resonant Feynman diagrams may significantly af- 
fect the u e u e HH cross section, and one may worry whether the averaging procedure employed 
here does yield credible results. To justify the averaging procedure, I compare in Fig. 10 
the SM e+e" -> v^^vtjjH, H -> WW -> 4j M 6j distribution for m H = 180 GeV and 
y/s = 3 TeV (black dashed line) with the result obtained from averaging the v e v e l^V(W H , 
W -> jj, H -> -> 4j (blue line) and u e u e WjjH, W ^ lv e H WW -> 4j (black 

solid line) cross sections (red line). For comparison, the magenta line shows the v e v e HH , 
if if — > 41T^ — > £^6j signal cross section. The red and black dashed lines agree within a few 
percent, lending credibility to the averaging procedure used in Fig. 9. 

Adopting the averaging procedure introduced above, I show in Fig. 11 the 6j invariant 
mass distribution for the SM e + e~ — > vvHH — > vvlvfoj signal (solid black line), AXhhh — 
+ 1 (dashed black line), and AX HHH = — 1 (dotted black line) for ^Js = 1 TeV and 3 TeV, 
together with the estimated uuiufij cross section including non-resonant diagrams (magenta 
lines). The non-resonant diagrams are seen to somewhat reduce the cross section over most 
of the Qj invariant mass range. Only for large values of M 6 j do they increase the differential 
cross section. The blue hatched histogram in Fig. 11a shows the W6j background. For 
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M 6j (GeV) 

FIG. 9. The SM e + e~ — * v e v e HH — * v e v e l~^vn§j cross section as a function of the 6j invariant 
mass, Mqj, for y/s = 1 TeV. The red line shows the e + e~ — > v e vjc .vijjH , H — > W + W~ — * Aj 
cross section. The blue curve gives the e + e~ — * v e D e AjH, H — > W + W~ — > Iv^jj cross section. The 
magenta line gives an estimate of the e + e~ — ► VeVel^vfij cross section, obtained from averaging 
the e+e" -» u^J^jjH, H -» VF+W" -> 4j and e+e" -> v e v e AjH, H -> -» 
differential cross sections. The cuts imposed and the efficiencies used are summarized in Eqs. (6), 
(8), (10), and (16). 

= 3 TeV the W6j cross section is too small to show up for the range of cross sections 
displayed. Clearly, the e + e~ — ► W6j background is negligible at both center of mass energies. 

If all four W bosons in HH — > AW decay hadronically, the final state consists of eight 
jets and missing transverse momentum. To identify vvHH events in the $ T 8j final state, I 
require, in addition to the standard j> T and jet identification cuts of Eq. (6), four jet pairs 
which satisfy Eq. (10). Furthermore, the jet pairs are required to form two groups of four 
jet systems which satisfy Eq. (16). The four jets have to consist of two jet pairs with each 
jj system fulfilling Eq. (10). 

The main background to the $ T 8j final state originates from non-resonant diagrams which 
I estimate by calculating the cross section for e + e~ — ► uuAjH with H — > WW — ► Aj and 
employing the same averaging procedure as for the £$ T 6j final state. The background from 
e + e~ — > 8j with the missing transverse momentum originating from jet mismeasurements 
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FIG. 10. The SM e + e~ — ► v e v e {~V£jjH, H — > WW — > Aj cross section as a function of the 
6j invariant mass, Mqj for y/s = 3 TeV (dashed black line). The black solid line shows the 
v e v e W jjH , W — > £vj>, H — > WW — > 4j cross section. The blue curve gives the v^^viWH, 
W — » jj, H — » VFW^ — ► 4j cross section. The black dashed line gives an estimate of the 
e + e _ — » v e v e t~^vijjH , — ► — > 4j cross section cross section from averaging the v e v e W jjH , 
W -+lv t ,H -> WW -> 4j and fe^i/^ff, W jj, if WW 4j differential cross sections. 
For comparison, the magenta line shows the v e v e HH, HH — > 4PF — > £f£6j signal cross section. 
The cuts imposed and the efficiencies used are summarized in Eqs. (6), (8), (10), and (16). 

is expected to be small except for Higgs pair invariant masses close to \fs. The results 
for the Higgs pair invariant mass distribution with ^/s = 1 TeV and 3 TeV are shown in 
Fig. 12. The non-resonant diagrams included in e + e _ — > uuAjH, H — > WW — > 4j, are seen 
to substantially reduce the cross section away from the threshold (magenta line). There 
is no guarantee that the averaging procedure used approximates the ui>8j cross section 
including the full set of non-resonant Feynman diagrams with a similar accuracy as that 
observed for the vviv^Sj final state. The non-resonant diagrams ignored in e + e~ — > uuAjH, 
H — > WW — > 4j may increase, or further decrease, the cross section. The main conclusion 
drawn from Fig. 12, therefore, is that non-resonant Feynman diagrams may substantially 
affect the uP8j cross section, and alter the shape of the Mjjh distribution. The effect of 
non-resonant Feynman diagrams can of course be reduced by imposing a more stringent cut 
on the 4j invariant mass than that used here (see Eq. (16)). Whether this will be possible 
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FIG. 11. The e + e — > vuHH — * vvtvfiy cross section as a function of the 6j invariant mass, 
M 6j , for m H = 180 GeV and a) y/a = 1 TeV and b) y/s = 3 TeV. The black solid line is the 
prediction of the SM signal cross section. The black dashed and dotted lines correspond to the 
uvHH signal cross section for AXhhh = +1 and AXhhh = — 1, respectively. The magenta 
line shows the estimated SM cross section for e + e~ — > vvtv^j. The blue hatched histogram in 
part a) shows the SM e + e~ — ► W6j background. The cuts imposed and the efficiencies used are 
summarized in Eqs. (6), (8), (10), and (16). 
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FIG. 12. The e + e~ — * uuHH — > vv$j cross section as a function of the Higgs pair invariant 
mass, for = 180 GeV and a) y/s = 1 TeV and b) y/s = 3 TeV. The black solid line is the 
prediction of the SM signal cross section. The black dashed and dotted lines correspond to the 
vvHH signal cross section for AXhhh = +1 and AXhhh = — 1, respectively. The magenta line 
shows the estimated SM cross section for e + e~ — > vv'&j obtained from the vvAjH, H — > WW — ► 4j 
cross section. The cuts imposed and the efficiencies used are summarized in Eqs. (6), (8), (10), 
and (16). 
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TABLE II. Cross sections in fb for e + e" — »■ vvbbbb for m H = 120 GeV and m H = 140 GeV with 
three or more 6-tags, e + e~ — > vvAjbb for m# = 140 GeV, and e + e~ — > vvtvfij and e + e~ — ► i/i/8j 
for rriH = 180 GeV. Results are shown for y/s = 1 TeV and 3 TeV, the SM signal, the full set of 
Feynman diagrams (labeled "all"), and the reducible backgrounds (labeled "bgd"). The reducible 
background does not include contributions where the missing transverse momentum originates only 
from the energy loss of 6-quarks and/or jet mismeasurements. The cuts and efficiencies used are 
listed in Eqs. (6), (8), (10), (11), (15), and (16). 



e + e — ► vvbbbb 





signal 


all 


bgd 




1 TeV, m# = 


120 GeV 


0.026 


0.024 


0.003 


V~s = 


1 TeV, rriff = 


140 GeV 


0.0044 


0.0052 


0.0014 


V~s = 


3 TeV, in h = 


120 GeV 


0.20 


0.16 


0.01 


V~s = 


3 TeV, 77i H = 


140 GeV 


0.040 


0.037 


0.006 


e + e — > vvAjbb 




signal 


all 


bgd 


V~s = 


1 TeV, mH = 


140 GeV 


0.0039 


0.0032 






3 TeV, rriff = 


140 GeV 


0.025 


0.023 




e + e — > vvivfij 




signal 


all 


bgd 


V~s = 


1 TeV, ttih = 


180 GeV 


0.0065 


0.0054 


4.0 x 10" 5 




3 TeV, ttih = 


180 GeV 


0.050 


0.041 


< 1 x 10~ 6 


e + e — > 1/1/8,7 




signal 


all 


bgd 


V~s = 


1 TeV, m H = 


180 GeV 


0.0047 


0.0036 






3 TeV, ttih = 


180 GeV 


0.032 


0.022 





largely depends on the resolution of the hadronic calorimeter which can be achieved in an 
ILC/CLIC experiment. 

D. Compilation of cross sections 

Before I derive sensitivity limits for the Higgs boson self-coupling, I present, in Table II, 
integrated signal and background cross sections for e + e~ — > vvHH and the Higgs boson 
masses and final states discussed in Sees. VA - VC. Table II shows that, in contrast to 
ZHH production, the background is always small in e + e~ — > vvHH . Furthermore, the 
vvHH — > vvbbbb cross section is considerably larger than that for ZHH — > jjbbbb for 
TTiu — 120 GeV and m# = 140 GeV. Non-resonant diagrams have a moderate effect except 
for vvHH — > vvbbbb production and ttih = 140 GeV, and for e + e~ — > vv8j. For vvbbbb 
production, QCD diagrams contribute at 0(a^a 4 ). Since these diagrams contribute only 
10 — 30% of the cross section over most of the M H h range, I estimate the renormalization 
scale uncertainty of the vvHH — > vvbbbb cross section to be no more than 10%. 

The main uncertainty of the e + e~ — > vvAjbb, e + e~ — > vvivfij and e + e~ — > vv8j cross 
sections originates from the unknown effect of non-resonant diagrams. Unfortunately, calcu- 
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lations of these processes including the full set of contributing Feynman diagrams is currently 
beyond what automated matrix element based programs can handle. I have presented re- 
sults based on approximations which include subsets of non-resonant diagrams, and argued 
that these should account for most of the effects of non-resonant diagrams. Nevertheless, 
an uncertainty of 0(20%) (0(50%)) remains for the e + e~ — > uuAjbb and e + e~ — > uuivfij 
(e + e~ — > vu8j) cross section. 



VI. SENSITIVITY LIMITS FOR AX HHH 

I now present quantitative sensitivity limits for the Higgs boson self-coupling for e + e~ — > 
ZHH and e + e~ — > vvHH , and the final states discussed in Sees. IV and V. Limits are 
derived from the M HH distribution except for e + e~ — > vufafij where the 6j invariant 
mass distribution is analyzed. Results are presented for integrated luminosities of 0.5 ab _1 , 
1 ab -1 , and 2 ab' 1 at y/s = 0.5 TeV and 1 TeV (ILC), and 1 ab _1 , 2 ab" 1 , and 3 ab" 1 
at \/s = 3 TeV (CLIC). According to the current ILC (CLIC) design parameters [32,33], 
an integrated luminosity of 1 ab -1 (3 ab -1 ) corresponds to 5 years of running. For ZHH 
production, I analyze the jjbbbb final state with > 3 and 4 6-tags. For y/s = 0.5 TeV, 
bounds are calculated for the sets of efficiencies listed in Eqs. (8) and (9), otherwise I only 
use the efficiencies of Eq. (8). To derive limits, I use the cross sections obtained including 
non-resonant Feynman diagrams. 

As the statistical tool of choice I adopt a log likelihood test. The expression for the 
log-likelihood function is 



-21ogL = -2 



E (-fsSi - f B B t + n 0l \og(f s S l + f B Bi) - log(n 0l !)) 



I (fs - I) 2 (fs - If ( , 

+ (Afsf + (A/ B ) 2 ' 1 ) 

The sum extends over the number of bins, Si and Bi are the number of signal and background 
events in the ith bin, and noi is the number of reference (eg. SM) events in the ith bin. The 
uncertainties on the signal and background normalizations are taken into account via two 
multiplicative factors, fs and which are allowed to vary but are constrained within the 
relative uncertainties of the signal and background cross sections, A/5 and A/b, respectively. 

In order to simplify the analysis, I assume a common uncertainty for signal and back- 
ground, fs = /b = /, in the following. This can be justified by noting that either the 
signal is considerably larger than the background, or vice versa. In the first case, I use the 
numerical value of fs, in the second, /b is used. The uncertainties, A/, are determined 
individually from the renormalization scale uncertainty (if QCD diagrams contribute) or the 
approximation used to calculate the cross section. These uncertainties have been discussed 
in Sees. IV and V. In all other cases, a generic theory uncertainty of 10% is assumed to 
account for unknown higher order electroweak corrections. The values of A/ used in the 
following analysis are collected in Table III. 

The rather large uncertainties listed for \fs = 0.5 TeV, m H = 140 GeV, and/or > 3 b- 
tags originate from the large renormalization uncertainty of the (reducible) background and 
could be reduced by a NLO QCD calculation of e + e~ — > bbcjjj and e + e~ — > bbccjj. Such 
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TABLE III. Theoretical uncertainties, Af, used in the statistical analysis. 



e + er -► jjbbbb, e b = 0.9, P c ^ b = 0.1, P^ b = 0.005 



m H = 120 GeV 


4 6-tags > 3 6-tags 


m fl = 140 GeV 


4 6-tags > 3 6-tags 


v/i = 0.5 TeV 
y/s = 1 TeV 


10% 100% 
10% 10% 


y/s = 0.5 TeV 
y/s = 1 TeV 


40% 100% 
10% 30% 


e+e" -► jj6666, e fe = 0.8, P c ^ b = 0.02, = 0.001 


m H = 120 GeV 


4 6-tags > 3 6-tags 


m H = 140 GeV 


4 6-tags > 3 6-tags 


V« = 0.5 TeV 


10% 40% 


y/s = 0.5 TeV 


10% 100% 


e + e — > vvbbbb 


m H 


y/s = 1 TeV v/s = 3 TeV 


m H 


y/s = l TeV y/s = 3 TeV 


120 GeV 


10% 10% 


140 GeV 


10% 10% 


e + e — > uvAjbb, mu = 140 GeV 


e + e — > vvtvfoj, mu = 180 GeV 


e + e~ — > vD8j, mu = 180 GeV 


>/s = 1 TeV ^ = 3 TeV 


v/s = 1 TeV ^/i = 3 TeV 


v/s = 1 TeV ^ = 3 TeV 


20% 20% 


20% 20% 


50% 50% 



calculations are beyond the current state of the art of one-loop calculations. Likewise, the 
uncertainties listed for e + e~ — * uu4jbb, e + e~ — > uuiufij and e + e~ — > z/z/8j could potentially 
be reduced by performing a full matrix element based tree level calculation of these processes. 
Further advances in automated tree level programs may make this possible. Alternatively, 
independent measurements of the relevant cross sections away from the signal region could 
be used to reduce the theoretical uncertainties. 

If /s = /b = /, logL can be minimized analytically and one finds the minimum of logL 
to occur at 

/ = l - (l - (AffN + V / (l-(A/)2iV)2 + 4(A/)2iV ) , (18) 

where 

N = J2(S l + B l ) (19) 

i 

is the total number of events, 

No = E"oi (20) 

i 

the total number of reference events, and Af is the uncertainty of the reference cross section. 

The 68.3% confidence level (CL) limits which can be achieved in e + e~ — > ZHH — > jjbbbb 
are listed in Table IV. For a Higgs boson with mass close to the current lower mass bound [1] , 
^hhh can be measured with a precision of 30 — 60% in e + e~~ — > jjbbbb at y/s = 0.5 TeV 
with an integrated luminosity of 0.5 — 2 ab~ x , if one requires 4 6-tags. This result is in 
qualitative agreement with that reported in Ref . [53] . Since the ZHH cross section falls with 
increasing center of mass energy, the sensitivities which can be achieved for the same final 
state and Higgs boson mass at y/s = 1 TeV are slightly worse. While Xhhh can be measured 
with reasonable precision in jjbbbb production for mu = 120 GeV and y/s = 0.5 TeV, the 
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TABLE IV. Sensitivities achievable at 68.3% CL for the Higgs boson self-coupling, AXhhh 
(see Eq. (13)), in e+e" -► ZHH -► jjbbbb for m H = 120 GeV and 140 GeV, yfs = 0.5 TeV and 
1 TeV, and several choices of integrated luminosities. Results are shown for 4 and > 3 6-tags and, 
for = 0.5 TeV, two choices of 6-tagging efficiencies and misidentification probabilities. I assume 
the uncertainties listed in Table III. The cuts imposed are described in Sees. Ill and IV. 



e b = 0.9, P c 


= 0.1, P^ 6 = 


0.005 






0.5 ab- 1 


1 ab" 1 


2 ab" 1 


yfe = 0.5 TeV, m H = 120 GeV, 4 6-tags 


+0.60 
-0.56 


+0.44 
-0.41 


+0.33 
-0.30 


y/s = 0.5 TeV, m H = 120 GeV, > 3 6-tags 


+1 4 
-5.2 


_i_n qq q i 
-1.12 ' -4.8 


_i_n 70 —3 5 
-0.75 ' -4.5 


y/s = 0.5 TeV, m H = 140 GeV, 4 6-tags 


+2.6 
-7.7 


+2 1 
-6.8 


+1.7 -3.3 
-2.0' -6.0 


y/s = 0.5 TeV, m H = 140 GeV, > 3 6-tags 


+4.8 
-8.9 


+3.6 
-7.7 


+2.7 
-6.9 


y/s = l TeV, m H = 120 GeV, 4 6-tags 


+0.76 
-0.60 


+0.53 
-0.45 


+0.38 
-0.33 


y/s = l TeV, m H = 120 GeV, > 3 6-tags 


+0.65 
-0.58 


+0.46 
-0.42 


+0.33 
-0.31 


y/s = l TeV, m H = 140 GeV, 4 6-tags 


+1.6 
-4.2 


+1.1 -2.4 
-1.0 ' -3.5 


+0.77 
-0.71 


y/s = l TeV, = 140 GeV, > 3 6-tags 


+1.4 
-4.1 


+1.1 -2.0 
-1.1 ' -3.5 


+0.78 
-0.77 


e b = 0.8, P c ^ b = 0.02, P^ 6 = 0.001, > 3 6-tags 




0.5 ab" 1 


1 ab" 1 


2 ab" 1 


y/s = 0.5 TeV, m H = 120 GeV 


+1.0 -4.2 
-0.8 ' -4.9 


+0.80 
-0.65 


+0.60 
-0.51 


v/i = 0.5 TeV, m H = 140 GeV 


+5.1 
-8.6 


+3.7 
-7.3 


+2.7 
-6.4 



significantly reduced signal cross section and the increased background make a measurement 
of the Higgs boson self-coupling very difficult for m# = 140 GeV in this channel. The 
sensitivities can be improved somewhat by including other final states such as HH — > 4jbb, 
however, the resulting limits are still extremely loose. 

The bounds presented here for m H = 140 GeV and y/s = 0.5 TeV are considerably worse 
than those derived from the cross section analysis presented in Ref . [26] . This can be traced 
to the large non-resonant background, and the considerable theoretical uncertainty in the 
cross section. Both were not taken into account in the earlier analysis. 

At a 1 TeV machine, the sensitivity limits which can be achieved in the jjbbbb final 
state with 4 6-tags and m H = 140 GeV are approximately a factor 2 less stringent than 
those found for m# = 120 GeV. Although the signal cross section increases by about a 
factor 1.5 including the final state with 3 tagged 6-quarks, the much increased reducible 
background, combined with the substantial renormalization uncertainty of the background, 
ruin the gain from the increased signal cross section for y/s = 0.5 TeV. For y/s = 1 TeV, 
a slight improvement in the sensitivity limits is observed by including the final state with 
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three tagged 6-quarks. 

A straightforward technique to reduce the (reducible) background is to choose a 6-tagging 
efficiency somewhat smaller than that of Eq. (8) for which the charm and light quark/gluon 
mistagging probabilities are lower. Table IV shows that choosing the parameters of Eq. (9) 
does indeed improve the limits which can be achieved for y/s = 0.5 TeV and > 3 fe-tags, 
however, the gain is not sufficient to compensate for the increased background and theoretical 
(renormalization scale) uncertainty which results from including final states with 3 6-tags. 
Nevertheless, a dedicated search for the 6-tagging efficiency and charm and light quark/gluon 
mistagging probability which optimizes the sensitivity limits for Xhhh m &y prove beneficial. 
For a first step in this direction, see Ref. [28]. For y/s = 1 TeV, the background is always 
relatively small and not much is gained by varying e& and the charm and light quark/gluon 
jet misidentification probabilities. 

One of the reasons for the weak sensitivity bounds for > 3 tagged 6-quarks is the large 
renormalization uncertainty of the background cross section. Reducing the renormaliza- 
tion uncertainty requires either the calculation of the NLO QCD corrections for the main 
background sources, bbcjjj, bb4j and bbccjj production, or a precise measurement of the 
cross section of the processes. It is interesting to investigate how much the sensitivity limits 
would actually improve if the normalization uncertainty of the background could be reduced 
to Af = 10%. The sensitivity bounds which one can hope to achieve for y/s = 0.5 TeV, > 3 
tagged 6-quarks and the uncertainties listed in Table III, and those obtained for Af = 10%, 
are compared in Table V. While a reduced theoretical uncertainty will substantially improve 
the sensitivity limits for > 3 tagged 6-quarks, more stringent limits than those found for 
4 tagged 6-quarks are only found for ran = 140 GeV and the efficiencies of Eq. (9). Of 
course, further improvements may be possible by optimizing the 6-tagging efficiency and 
charm and light quark/gluon misidentification probabilities. The sensitivity limits which 
can be achieved for Xhhh i n e + e~ — > vvHH with the theoretical uncertainties of Table III 
are listed in Table VI. For van = 140 GeV (m# = 180 GeV), the combined limits from 
e + e~ — > vvbbbb and e + e~ — > uuAjbb (e + e~ — > vvivfij and e + e _ — > vv&j) are shown. Since 
the bbbb and bbjj backgrounds do not affect the Mhh differential cross section in the region 
sensitive to the Higgs boson self-coupling (see Sees. VA and VB), they have not been taken 
into account in the analysis. The bounds on the Higgs self-coupling which can be achieved 
in e + e~ — > vvHH for itlh = 120 GeV and 140 GeV at \/s = 1 TeV are seen to be consid- 
erably more stringent than those found for e + e _ — > ZHH (see Table IV). For ^/s = 3 TeV, 
Xhhh can be measured with a precision of 10 — 20% for the Higgs boson masses and the 
integrated luminosities considered here. If the theoretical uncertainty on the cross section 
for m H = 180 GeV can be reduced to 10%, the bounds listed in Table VI for y/s = 3 TeV 
improve by a factor 1.6 — 2. For itlh = 140 GeV (run = 180 GeV), the sensitivity limits are 
expected to improve by roughly a factor 1.1 (1.2) if additional final states are included in 
the analysis. 

The bounds derived in this Section should be compared with those one hopes to achieve 
at the LHC, a luminosity upgraded LHC (SLHC) and a Very Large Hadron Collider 
(VLHC) [19, 20] 2 . For the Higgs mass range allowed by current experimental data [1-3], 



2 Ref. [19,20] quotes 95% CL limits. I have recalculated the sensitivity limits for 68.3% CL with 
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TABLE V. Sensitivities achievable at 68.3% CL for the Higgs boson self-coupling, AXhhh 
(see Eq. (13)), in e+e" -> ZHH -> with > 3 6-tags for m H = 120 GeV and 140 GeV, 

yfs = 0.5 TeV, and several choices of integrated luminosities. Results are shown for the uncer- 
tainties listed in Table III and for Af = 10%, and for two choices of 6-tagging efficiencies and 
misidentification probabilities. The cuts imposed are described in Sees. Ill and IV. 



e b = 0.9, P c ^ b = 0.1, Pj^ b = 0.005 





0.5 ah" 1 1 ah" 1 2 ab" 1 


m H = 120 GeV, Af = 100% 
m H = 120 GeV, Af = 10% 


+1.4 +0.99 -3.1 +0.70 -3.5 
-5.2 -1.12 ' -4.8 -0.75 ' -4.5 

+0.86 -4.1 +0.66 -4.2 +0.52 
-1.04 ' -5.1 -0.74 ' -4.7 -0.56 


m H = 140 GeV, Af = 100% 
m H = 140 GeV, A/ = 10% 


+4.8 +3.6 +2.7 
-8.9 -7.7 -6.9 

+3.5 +2.8 +2.3 
-8.0 -7.5 -6.8 


e b = 0.8, P c ^ b = 0.02, Pj^ 6 = 0.001 




0.5 ah" 1 1 ab" 1 2 ab" 1 


m H = 120 GeV, A/ = 40% 


+1.0 -4.2 +0.80 +0.60 
-0.8 ' -4.9 -0.65 -0.51 


= 120 GeV, A/ = 10% 


+0.63 +0.47 +0.36 
-0.63 -0.46 -0.35 


m H = 140 GeV, Af = 100% 


+5.1 +3.7 +2.7 
-8.6 -7.3 -6.4 


m H = 140 GeV, A/ = 10% 


+2.5 +1.9 +1.5 
-7.7 -6.9 -6.3 


the LHC will not be able to probe the Higgs boson self-coupling. At the SLHC, with an 



integrated luminosity of 6 ab -1 , Xhhh can be determined with an accuracy of 50 — 70% 
for ma = 120 GeV from HH — > 6677, and 10 — 15% for ran = 180 GeV using the 
HH — >■ 414^ — > £ ± £' ± ^ T 4j final state. There are not enough HH — > 6677 events for 
m H = 140 GeV for a viable statistical analysis at the SLHC. A luminosity upgrade of the 
LHC may be realized at roughly the same time or earlier than an ILC with y/s = 0.5 TeV. 
Tables III and IV show that, for m# = 120 GeV, the ILC can measure the Higgs boson 
self-coupling with considerably higher precision than the SLHC. A measurement of Xhhh 
for m H = 140 GeV will be very difficult at both machines. For m H = 180 GeV, the SLHC 
shows a clear advantage, however, a recent combination of electroweak precision data and 
direct limits from the Tevatron experiments excludes a SM Higgs boson with m H = 180 GeV 
at more than 2a [3]. 

An upgraded ILC with \fs = 1 TeV considerably improves the chances to determine 
^hhh for m H = 140 GeV (see Tables IV and VI). For m H = 180 GeV, however, the 
SLHC still promises a more precise measurement. Since sensitivity bounds scale roughly 
like (/ Cdt)^ 4 \ this statement is true as long as an integrated luminosity larger than 2 ab -1 



otherwise unchanged parameters. The 68.3% CL limits are used for the comparison presented here. 
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TABLE VI. Sensitivities achievable at 68.3% CL for the Higgs boson self-coupling, AXhhh 
(see Eq. (13)), in e + e~~ — > vvHH for ^/s = 1 TeV and 3 TeV and several choices of integrated 
luminosities. For tuh = 120 GeV, results are shown for the vvbbbb final state with > 3 tagged 
6-quarks. For tuh = 140 GeV, the combined limits from e + e~ — > vvbbbb with > 3 6-tags and 
e + e~ — > vvAjbb are listed, and for m# = 180 GeV the combined limits from e + e~ — > vvlvfij 
and e + e~ — > z^8j are given. The cuts imposed are described in Sees. Ill and V. The theoretical 
uncertainties used are listed in Table III. 



y/s = 1 TeV 





0.5 ah" 1 


1 ab^ 1 


2 ab^ 1 


m H = 120 GeV 


+0.58 
-0.27 


+0.30 
-0.21 


+0.20 
-0.17 


m H = 140 GeV 


+0.99 
-0.41 


+0.94 
-0.38 


+0.78 
-0.25 


m H = 180 GeV 


+0.56 
-0.32 


+0.55 
-0.29 


+0.59 
-0.28 


Vs = 3 TeV 




1 ab" 1 


2 ab" 1 


3 ab" 1 


m H = 120 GeV 


+0.14 
-0.12 


+0.11 
-0.10 


+0.10 
-0.09 


m H = 140 GeV 


+0.15 
-0.19 


+0.15 
-0.15 


+0.11 
-0.14 


m H = 180 GeV 


+0.16 
-0.20 


+0.15 
-0.13 


+0.12 
-0.12 



can be accumulated at the SLHC. 

The results for a 3 TeV e + e~ collider (CLIC) should be compared with those which 
one can hope to achieve at a VLHC. Such a machine, assuming a center of mass energy of 
y/s = 200 TeV and an integrated luminosity of 600 fb -1 , can determine the Higgs boson 
self-coupling for itlh = 120 GeV (mn = 140 GeV) in the HH — > bbyy final state with a 
precision of 30 — 50% (40 — 60%). At CLIC, on the other hand, a precision of 10 — 20% can 
be achieved for both values of m H . For m H = 180 GeV, a VLHC can achieve a precision of 
about 2%. Assuming that a 10% theoretical uncertainty on the cross section can eventually 
be realized, the sensitivity of a 3 TeV e + e~ collider obtained from the combined vvlv^Sj and 
uu8j final states is limited to 7 — 12%. The two final states together account for about 43% 
of all HH — > AW decays. Including more 4W final states thus is expected to improve these 
limits by roughly a factor 2 1 / 4 1.2. A VLHC therefore may yield more precise information 
on the Higgs boson self-coupling than CLIC for m H = 180 GeV. 

VII. SUMMARY AND CONCLUSIONS 

After discovery of an elementary Higgs boson at the LHC, and tests of its fermionic and 
gauge boson couplings, experimental evidence that the shape of the Higgs potential has the 
form required for electroweak symmetry breaking will complete the proof that both, fermion 
and weak boson masses, are generated by spontaneous symmetry breaking. To probe the 
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Higgs potential, one must determine the Higgs boson self-coupling. 

Only Higgs boson pair production at colliders can accomplish this. Several years ago, 
studies have appeared in the literature [18-22], exploring the potential of the LHC, a lu- 
minosity upgraded LHC, and a Very Large Hadron Collider to probe the Higgs boson self- 
coupling. There are also numerous analyses of Higgs boson pair production at e + e~ collid- 
ers [16,17,23-30] which usually only explore one particular Higgs mass and/or final state, or 
only one center of mass energy. Furthermore, backgrounds are estimated using a leading-log 
shower approximation, and the effect on non- resonant diagrams is not taken into account. 

In this paper, I have presented a more comprehensive analysis of Higgs pair production 
in e + e~ collisions. Both ZHH and vvHH production have been investigated for several 
Higgs boson masses and center of mass energies. The cross section for vvHH production 
is much smaller than that for ZHH production for y/s = 0.5 TeV and below. However, it 
grows quickly with increasing energies, and becomes the dominant source of Higgs boson 
pairs for y/s > 1 TeV. ZHH (vuHH) production therefore was studied for y/s < 1 TeV 
{\fs > 1 TeV) only. Acceptance cuts and minimal detector effects in form of Gaussian 
smearing, as well as the energy loss of 6-quarks were taken into account in all calculations. 
The efficiencies and misidentification probabilities used are summarized in Eqs. (8) and (9). 
The distribution of the Higgs pair invariant mass, Mhh, was found to be sensitive to the 
Higgs boson self-couplings and was used for a log-likelihood based sensitivity analysis. 

For ZHH production, I concentrated on the ZHH — > jjbbbb final state, requiring that 
the two light quark jets are consistent with a Z-boson. Requiring four tagged 6-quarks, 
the jjbbcc and bb4j backgrounds were found to be relatively small for itlh = 120 GeV and 
y/s = 0.5 TeV. Non-resonant diagrams change the jjbbbb cross section by about 10%. For 
mu — 140 GeV, on the other hand, taking into account the non- resonant diagrams roughly 
doubles the cross section, while the jjbbcc background can easily be as large as the signal 
cross section. This, together with the extremely small signal cross section, and the large 
renormalization uncertainty of the reducible background cross section, implies that it will 
be extremely difficult to measure the Higgs boson self-coupling in jjbbbb production with 
four tagged 6-quarks and ran = 140 GeV at an ILC with a center of mass energy of 0.5 TeV. 
At y/s = 1 TeV, the reducible and non-resonant backgrounds are significantly smaller than 
at a 0.5 TeV collider, and a rough measurement of Xhhh in e + e~ — > jjbbbb may be possible 
for m H = 140 GeV. 

To increase the signal cross section, one can include other Z and/or Higgs decay final 
states. Since BiZ — > £ + £~) is small, there will be little gain by including the l + l~bbbb final 
state in the analysis. For m H = 140 GeV, B(H -> WW* -> 4j) « 23%. If one of the Higgs 
bosons in ZHH production decays into four jets, the final state consists of bb6j. Since the 
combinatorial background complicates identification of which jets originate from Higgs and 
which from Z-boson decays, I have not considered bb6j production here. 

A more straightforward approach is to increase the signal cross section by taking into 
account final states where only three 6-quarks are tagged. Requiring only three tagged b- 
quarks strongly increases the reducible background and this more than compensates the 
potential gain in sensitivity from the increased signal cross section. However, the reducible 
background can be reduced by choosing a different working point in — P c ^b — Pj^b space. 
Replacing the efficiencies and misidentification probabilities of Eq. (8) with those of Eq. (9), 
however, yields only a minor improvement in the sensitivity to A. Nevertheless, a search for 
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an optimal working point similar to that in Ref. [28] may be promising. 

For y^s = 1 TeV, the vvbbbb cross section is about a factor two larger than that for jjbbbb 
production. Furthermore, the reducible background is small, both for mu = 120 GeV and 
140 GeV. Non-resonant diagrams modify the M HH differential cross section by 0(10%) 
except close to the HH threshold where they have a much larger effect. The effect from 
non-resonant diagrams therefore can be easily confused with a positive anomalous Higgs 
self-coupling, Xhhh > (see Figs. 6a and 7a). Due to the larger cross section, and the 
reduced background, the bounds which can be obtained for m# = 120 GeV are significantly 
better than those one may hope to achieve in ZHH — > jjbbbb. 

For m H = 140 GeV, the signal cross section can be increased by including the vvbbAj final 
state which results when one of the Higgs bosons decays via H — > WW* — > Aj. Combining 
the limits from vvbbbb and vvbbAj production, I found that the Higgs boson self-coupling can 
be probed with 25 — 80% accuracy with 2 ab _1 at a 1 TeV e + e~ collider. Such a machine will 
also be able to probe A for a heavier Higgs boson. For m# = 180 GeV, and utilizing the two 
final states with the largest branching ratios, HH — > AW — > fi^vfij and HH — > AW — > 8j, 
I found that a 1 TeV e + e~ collider can probe A with an accuracy of 30 — 60%. A 3 TeV 
collider will be able to measure the Higgs boson self-coupling with a precision of 10 — 20% 
for the Higgs boson masses investigated here. 

A MadEvent based calculation of the full set of non-resonant Feynman diagrams con- 
tributing to vvbAj, vvtv&j and vv8j production is currently not feasible as it requires very 
large computing resources. Other programs, which are not based on the evaluation of Feyn- 
man diagrams promise a solution for this problem in the future, however, these programs 
are currently not able to handle the final states of interest. In absence of a calculation which 
includes the full set of non-resonant Feynman diagrams, I have estimated their effect from 
calculations which include sub-sets of the non-resonant Feynman diagrams. Unavoidably, 
this introduces a theoretical uncertainty in the cross section which I have assumed to be of 
the same size as the effect from those non-resonant diagrams which were included in the cal- 
culation. If this uncertainty can be reduced in the future by a calculation which includes the 
full set of non-resonant Feynman diagrams, the bounds which can be achieved will improve 
by a factor 1.6 - 2. 

Uncertainties in the calculations presented in this paper also originate from the detector 
performance assumed. A better resolution of the hadronic calorimeter than that assumed 
such as that promised by the particle flow algorithm [64] may make it possible to tighten 
the invariant mass window for 6-jet and light quark/gluon jet pairs. This would reduce 
both reducible and irreducible backgrounds and thus improve the sensitivity to anomalous 
Higgs boson self-couplings. The reducible backgrounds in final states involving 6-quarks 
strongly depend on the charm and light quark/gluon jet misidentification probabilities. More 
powerful algorithms for tagging 6-quarks may well reduce the misidentification probabilities 
in the future. 

One of the main results of the calculations presented here is that a 0.5 TeV ILC will be 
able to determine the Higgs boson self-coupling only if the Higgs boson mass is rather close 
to the current lower bound. In contrast, at a 1 TeV e + e~ linear collider it will be possible 
to measure A for larger values of run, and to reach a significantly higher precision than for 
yt~s = 0.5 TeV. The precision which can be reached at a 1 TeV ILC for < 140 GeV 
is considerably better than that one can hope to achieve at the SLHC or a VLHC. For a 
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heavier Higgs boson with m# = 180 GeV, which is disfavored by a global fit to electroweak 
data and current direct limits, on the other hand, the SLHC and VLHC promise better 
limits. Even a 3 TeV e + e~ collider will not be able to measure the Higgs boson self-coupling 
for m H = 180 GeV more accurately than a VLHC. 



ACKNOWLEDGMENTS 



I would like to thank J. Reuter for his help with WHIZARD and for useful discussions, and 
F. Krauss for an in-depth discussion on the capabilities of Sherpa. I would also like to thank 
the High Energy Physics Group at Michigan State University, where most of this work was 
done, for their generous hospitality. This research was supported by the National Science 
Foundation under grants No. PHY-0456681 and PHY-0757691. 



35 



REFERENCES 



[1] [ALEPH Collaboration and CDF Collaboration and DO Collaboration and DELPHI 
Collaboration and L3 Collaboration and OPAL Collaboration and SLD Collaboration 
and LEP Electroweak Working Group and Tevatron Electroweak Working Group and 
SLD Electroweak Working Group and Heavy Flavour Group ], arXiv:081 1.4682 [hep-ex]; 
H. Flaecher, M. Goebel, J. Haller, A. Hoecker, K. Moenig and J. Stelzer, Eur. Phys. J. 
C 60, 543 (2009) [arXiv:0811.0009 [hep-ph]]. 

[2] The TEVNPH Working Group, arXiv:0903.4001 [hep-ex]. 

[3] M. Goebel, arXiv:0905.2488 [hep-ph]. 

[4] M. Dittmar and H. K. Dreiner, Phys. Rev. D 55, 167 (1997); D. Rainwater and D. Zep- 
penfeld, Phys. Rev. D 60, 113004 (1999) [Erratum-ibid. D 61, 099901 (2000)]; N. Kauer, 
T. Plehn, D. Rainwater and D. Zeppenfeld, Phys. Lett. B 503, 113 (2001); N. Akchurin 
et al, CMS-NOTE-2002/066; B. Mellado, ATL-CONF-2002-004. 

[5] S. Asai et al, ATL-PHYS-2003-005; G. Azuelos and R. Mazini, ATL-PHYS-2003-004; 
M. Diihrssen, ATL-PHYS-2003-010. 

[6] ATLAS TDR, report CERN/LHCC/99-15 (1999); G. Aad et al. [The ATLAS Collabo- 
ration], arXiv:0901.0512. 

[7] G. L. Bayatian et al. [CMS Collaboration], J. Phys. G 34, 995 (2007). 

[8] D. Rainwater, D. Zeppenfeld and K. Hagiwara, Phys. Rev. D 59, 014037 (1999); 
T. Plehn, D. Rainwater and D. Zeppenfeld, Phys. Lett. B 454, 297 (1999) and Phys. 
Rev. D 61, 093005 (2000). 

[9] D. Zeppenfeld, R. Kinnunen, A. Nikitenko and E. Richter-Was, Phys. Rev. D 62, 
013009 (2000); M. Hohlfeld, ATL-PHYS-2001-004; M. Diihrssen, S. Heinemeyer, H. Lo- 
gan, D. Rainwater, G. Weiglein and D. Zeppenfeld, Phys. Rev. D 70, 113009 (2004) 
[arXiv:hep-ph/0406323]; R. Lafaye, T. Plehn, M. Rauch, D. Zerwas and M. Diihrssen, 
arXiv:0904.3866 [hep-ph]. 
[10] D. Rainwater, Phys. Lett. B 503, 320 (2001); V. Drollinger, T. Miiller and D. Denegri, 
arXiv:hep-ph/0201249. 

[11] V. Drollinger, T. Miiller and D. Denegri, arXiv:hep-ph/0111312; V. Kostioukhine, 
J. Leveque, A. Rozanov, and J.B. de Vivie, ATL-PHYS-2002-019; D. Green et al, 
FERMILAB-FN-705 (August 2001); F. Maltoni, D. Rainwater and S. Willenbrock, 
Phys. Rev. D 66, 034022 (2002); A. Belyaev and L. Reina, JHEP 0208, 041 (2002); 
A. Belyaev, F. Maltoni and L. Reina, in Proc. of the APS/DPF/DPB Summer Study on 
the Future of Particle Physics (Snowmass 2001) ed. N. Graf, arXiv:hep-ph/0110274. 

[12] O. J. Eboli and D. Zeppenfeld, Phys. Lett. B 495, 147 (2000). 

[13] T. Plehn and D. Rainwater, Phys. Lett. B 520, 108 (2001); T. Han and B. McElrath, 
Phys. Lett. B 528, 81 (2002). 

[14] J. A. Aguilar-Saavedra et al. [ECFA/DESY LC Physics Working Group Collabora- 
tion], arXiv:hep-ph/0106315 and references therein; T. Abe et al. [American Linear 
Collider Working Group Collaboration], in Proc. of the APS/DPF/DPB Summer Study 
on the Future of Particle Physics (Snowmass 2001) ed. R. Davidson and C. Quigg, 
arXiv:hep-ex/0 106056 and references therein; A. Djouadi et al. [ILC Collaboration], 
arXiv: 0709. 1893 [hep-ph]. 

[15] D. A. Dicus, C. Kao and S. S. Willenbrock, Phys. Lett. B 203, 457 (1988); 
E. W. N. Glover and J. J. van der Bij, Nucl. Phys. B 309, 282 (1988); E. W. N. Glover 



36 



and J. J. van der Bij, CERN-TH-5022-88, in Proceedings of the "23rd Rencontres 
de Moriond: Current Issues in Hadron Physics, Les Arcs, France, Mar 13-19, 1988; 
G. Cynolter, E. Lendvai and G. Pocsik, hep-ph/0003008, Acta Phys. Polon. B 31, 1749 
(2000). 

[16] F. Boudjema and E. Chopin, Z. Phys. C 73, 85 (1996). 
[17] V. A. Ilyin et al, Phys. Rev. D 54, 6717 (1996). 

[18] F. Gianotti et al, Eur. Phys. J. C 39, 293 (2005) [arXiv:hep-ph/0204087]. 

[19] U. Baur, T. Plehn and D. Rainwater, Phys. Rev. Lett. 89, 151801 (2002) and Phys. 

Rev. D 67, 033003 (2003). 
[20] U. Baur, T. Plehn and D. Rainwater, Phys. Rev. D 69, 053004 (2004) [arXiv:hep- 

ph/0310056]. 

[21] U. Baur, T. Plehn and D. L. Rainwater, Phys. Rev. D 68, 033001 (2003) [arXiv:hep- 
ph/0304015]. 

[22] A. Blondel, A. Clark and F. Mazzucato, ATL-PHYS-2002-029 (November 2002). 

[23] A. Djouadi, W. Kilian, M. Miihlleitner and P. M. Zerwas, Eur. Phys. J. C 10, 27 (1999). 

[24] D. J. Miller and S. Moretti, Eur. Phys. J. C 13, 459 (2000). 

[25] M. Battaglia, E. Boos and W. M. Yao, in Proc. of the APS/DPF/DPB Summer Study 
on the Future of Particle Physics (Snowmass 2001) ed. R. Davidson and C. Quigg, 
arXiv:hep-ph/0111276. 

[26] C. Castanier, P. Gay, P. Lutz and J. Orloff, arXiv:hep-ex/0101028. 

[27] Y. Yasui et al, arXiv:hep-ph/0211047. 

[28] D. E. Boumediene and P. Gay, In the Proceedings of 2007 International Linear Collider 
Workshop (LCWS07 and ILC07), Hamburg, Germany, 30 May - 3 Jun 2007, pp HIG06 
[arXiv:0801.0866 [hep-ex]]. 

[29] Y. Takubo, arXiv:0901.3598 [hep-ph]. 

[30] M. F. Giannelli, arXiv:0901.4895 [hep-ex]. 

[31] R. Belusevic and G. Jikia, Phys. Rev. D 70, 073017 (2004) [arXiv:hep-ph/0403303]. 

[32] J. Brau et al, "International Linear Collider reference design report. 1: Executive sum- 
mary. 2: Physics at the ILC. 3: Accelerator. 4: Detectors," SLAC-R-857. 

[33] H. Braun et al. [CLIC Study Team Collaboration], CERN-OPEN-2008-021, CLIC- 
NOTE-764. 

[34] J. J. Lopez- Villarejo and J. A. M. Vermaseren, arXiv:0812.3750 [hep-ph]. 
[35] A. Gutierrez- Rodriguez, M. A. Hernandez- Ruiz, O. A. Sampayo, A. Chubykalo and 
A. Espinoza-Garrido, J. Phys. Soc. Jap. 77, 094101 (2008) [arXiv:0807.0663 [hep-ph]]. 
[36] G. Ambrosio et al, Report No. Fermilab-TM-2149 (2001). 

[37] B. W. Lee, C. Quigg and H. B. Thacker, Phys. Rev. Lett. 38, 883 (1977) and Phys. 

Rev. D 16, 1519 (1977). 
[38] T. Plehn and M. Rauch, Phys. Rev. D 72, 053008 (2005) [arXiv:hep-ph/0507321]. 
[39] S. Kanemura et al, Phys. Lett. B 558, 157 (2003). 

[40] W. Hollik and S. Penaranda, Eur. Phys. J. C 23, 163 (2002); A. Dobado, M. J. Herrero, 

W. Hollik and S. Penaranda, Phys. Rev. D 66, 095016 (2002). 
[41] D. B. Kaplan and H. Georgi, Phys. Lett. 136B, 183 (1984); H. Georgi, Phys. Lett. 

151B, 57 (1985). 

[42] C. Csaki et al., Phys. Rev. D 68, 035009 (2003); T. Han, H. E. Logan, B. McElrath and 
L. T. Wang, Phys. Rev. D 67, 095004 (2003); C. Dib, R. Rosenfeld and A. Zerwekh, 



37 



arXiv:hep-ph/0302068. 
[43] V. Barger et al. Phys. Rev. D 67, 115001 (2003). 

[44] F. Maltoni and T. Stelzer, JHEP 0302, 027 (2003) [arXiv:hep-ph/0208156]; J. Alwall 

et al, JHEP 0709, 028 (2007) [arXiv:0706.2334 [hep-ph]]. 
[45] M. L. Mangano, M. Moretti, F. Piccinini, R. Pittau and A. D. Polosa, JHEP 0307, 001 

(2003) [arXiv:hep-ph/0206293]. 
[46] T. Behnke et al. [ILC Collaboration], arXiv:0712.2356 [physics. ins-det]. 
[47] T. Behnke, S. Bertolucci, R. D. Heuer and R. Settles, DESY-01-011, March 2001. 
[48] S. M. Xella-Hansen, M. Wing, D. J. Jackson, N. de Groot, and C.J.S. Damerell, LC- 

PHSM-2003-061. 

[49] U. Baur and L. H. Orr, Phys. Rev. D 76, 094012 (2007) [arXiv:0707.2066 [hep-ph]]. 
[50] S. Hillert [LCFI Collaboration], In the Proceedings of 2005 International Linear Collider 

Workshop (LCWS 2005), Stanford, California, 18-22 Mar 2005, pp 0313. 
[51] G. Belanger et al, Phys. Lett. B 576, 152 (2003) [arXiv:hep-ph/0309010]; R. Y. Zhang, 

W. G. Ma, H. Chen, Y. B. Sun and H. S. Hou, Phys. Lett. B 578, 349 (2004) [arXiv:hep- 

ph/0308203]. 

[52] F. Boudjema et al, In the Proceedings of 2005 International Linear Collider Workshop 
(LCWS 2005), Stanford, California, 18-22 Mar 2005, pp 0601 [arXiv:hep-ph/0510184]. 

[53] T. Barklow, talk given at the 2007 International Linear Collider Workshop (LCWS 
2007), DESY, Hamburg, Germany, May 30 - June 3, 2007. 

[54] T. Barklow, talk given at the 2007 ALCPG Workshop, Fermilab, 22-26 October 2007. 

[55] T. Gleisberg et al. JHEP 0902, 007 (2009) [arXiv:081 1.4622 [hep-ph]]. 

[56] F. Krauss, R. Kuhn and G. Soff, JHEP 0202, 044 (2002) [arXiv:hep-ph/0109036]. 

[57] F. Krauss, private communication 

[58] T. Gleisberg and S. Hdche, JHEP 0812, 039 (2008) [arXiv:0808.3674 [hep-ph]]. 

[59] C. Duhr, S. Hoche and F. Maltoni, JHEP 0608, 062 (2006) [arXiv:hep-ph/0607057]. 

[60] W. Kilian, T. Ohl and J. Reuter, arXiv:0708.4233 [hep-ph]. 

[61] A. Cafarella, C. G. Papadopoulos and M. Worek, arXiv:0710.2427 [hep-ph]. 

[62] M. Moretti, T. Ohl and J. Reuter, arXiv:hep-ph/0102195. 

[63] K. Kolodziej and S. Szczypinski, Nucl. Phys. B 801, 153 (2008) [arXiv:0803.0887 [hep- 
ph]]; K. Kolodziej, arXiv:0903.3334 [hep-ph]. 
[64] S. R. Magill, New J. Phys. 9, 409 (2007). 



38 



